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ABSTRACT 

We investigate the possibility of probing dark energy by measuring the isotropy of 
the galaxy cluster autocorrelation function (an Alcock-Paczyhski test). The correla- 
tion function is distorted in redshift space because of the cluster peculiar velocities, 
but if these are known and can be subtracted, the correlation function measurement 
becomes in principle a pure test of cosmic geometry. Galaxy cluster peculiar velocities 
can be measured using the kinetic Sunyaev Zel'dovich (kSZ) effect. Upcoming CMB 
surveys, e.g., ACT, SPT, Planck, are expected to do this with varying levels of accu- 
racy, dependent on systematic errors due to cluster temperature measurements, radio 
point sources, and the primary CMB anisotropy. We use the Hubble volume iV-body 
simulation and the hydrodynamic simulation results of Nagai et. al (2003) to simulate 
various kSZ surveys. We find by model fitting that a measurement of the correlation 
function distortion can be used to recover the cosmological parameters that have been 
used to generate the simulation. However, the low space density of galaxy clusters 
requires larger surveys than are taking place at present to place tight constraints on 
cosmology. For example, with the SPT and ACT surveys, f^A could be measured to 
within 0.1 and 0.2 respectively at one sigma, but only upper limits on the equation of 
state parameter w will be possible. Nevertheless, with accurate measurements of the 
kSZ effect, this test can eventually be used to probe the dark energy equation of state 
and its evolution with redshift, with different systematic errors than other methods. 
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1 INTRODUCTION 

The Alcock-Paczyriski (1979, hereafter AP) test is a con- 
ceptually simple probe of cosmic geometry and hence dark 
energy. It is based on the fact that an intrinsically isotropic 
object will appear distorted when we translate it from ob- 
served space (angles and redshifts) into real space (units 
of length) if we use the wrong cosmology (to relate angles 
and redshift to sizes). The strength of the distortion will 
depend on cosmic geometry, for example the value of the 
cosmological constant A. The test is a way to measure dark 
energy which requires only that we find and measure known 
perfectly symmetric objects in the Universe. In the present 
work we use the autocorrelation function of galaxy clusters 
as the "isotropic" object. In addition to the possible geomet- 
ric distortions, the correlation function will have a distortion 
of different kind, which is due to peculiar velocities of clus- 
ters. This affects the line of sight component of cluster pair 
separations (e.g., as applied to galaxies by Davis and Pee- 
bles 1983, Kaiser 1987, and clusters by the Padilla & Baugh 
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2002). There is a way to remove this distortion however, 
by making use of the kinetic Sunyaev-Zeldovich (kSZ) effect 
(Sunyaev & Zeldovich 1972). The kSZ effect is the Doppler 
shift of CMB photons due to a hot intracluster medium. 
The CMB photons interact with free electrons and become 
Doppler shifted when the cluster is moving with respect to 
the CMB rest frame. This causes temperature fluctuations 
that depend on the bulk velocity of cluster. There are cur- 
rently on-going and future surveys that plan to measure this 
effect (see e.g., Vale and White 2005 for strategies). At the 
time of writing, the kSZ effect has not yet been detected for 
single clusters, although upper limits on have been estimated 
on the bulk velocity averaged over several clusters (e.g., 1420 
kms -1 at 95 % confidence by Benson et al. 2003). In the 
present paper we investigate the correction of redshift dis- 
tortions from peculiar velocities using kSZ measurements, 
and how this can leave the pure geometric distortions from 
which we can find cosmological parameters. 

Similar methods have been put forward for constraining 
cosmology with the clustering of different objects. Phillipps 
(1994) suggested using quasar clustering for the Alcock 
Paczyhski test (averaging over many close pairs of quasars) . 
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Several other studies have adopted the quasar correlation 
function as the isotropic object (e.g., Popowski et al. 1998, 
Hoyle et al. 2002, da Angela et al. 2005). The difficulty of 
this method is that we only observe the positions of quasars 
in redshift space, where the correlation function is distorted 
by peculiar velocities (e.g., Kaiser 1987, Hamilton 1992) and 
by other source of noise, such redshift errors. 

Many studies have suggested using the kinetic Sunyaev- 
Zel'dovich effect (e.g., Haehnelt & Tegmark 1996, Lange 
et al. 1998, Kashilinsky & Atrio-Barandela 2000, Aghanim 
et al. 2001, Atrio-Barandela et al. 2004, etc.) to measure 
peculiar velocities of galaxy clusters. There are several up- 
coming sky surveys that aim to do this, such as ACT 1 , SPT 2 
and Planck 3 In principle, these types of measurements can 
be used to eliminate the distortion due to peculiar veloci- 
ties and therefore we can recover the isotropic correlation 
function. The bulk cluster peculiar velocity is not the sole 
factor in the distortion and there are sources of noise in kSZ 
measurements. These include microwave background fluc- 
tuations, cluster internal velocities, and point source con- 
tamination (see simulations by Haehnelt & Tegmark 1996, 
Aghanim et al. 2005.) 

In the present work, we use large N-body simulations 
(the Hubble volume, Evrard et al. 2002), and the hydrody- 
namic simulation results of Nagai et al. (2003) to test how 
well the peculiar velocity-corrected galaxy cluster AP test 
can be carried out. Our plan for the paper is as follows: 
In §2 we discuss the theoretical cluster correlation function 
we use as a model as well as that measured from the Hub- 
ble volume simulation. In §3.1, we explain the model fitting 
procedure and predict the effect of kSZ measurements. We 
describe and measure the geometric distortions as a function 
of redshift. Using the recovered values we constrain cosmo- 
logical parameters from the simulation data, discussing the 
results in the same section. We show results from simula- 
tions of particular surveys (SPT and ACT) in §4.1-§4.2. We 
summarize our results and conclude in 55. 



2 CLUSTER CORRELATION FUNCTION 

The cluster autocorrelation function £(r) is a measure of the 
probability of finding a cluster at a distance r from another 
cluster; typically it is computed as the ratio of the number 
of clusters at a distance r from another cluster [DD(r)], 
divided by the number expected in the absence of clus- 
tering [DR(r)] minus 1, i.e. £(r) = — 1, where DD 
represent cluster-cluster pairs and DR cluster-random pairs. 
£(<7, 7r) carries the same definition but the distance between 
two clusters is represented in terms of the component di- 
rections perpendicular (a) and parallel to the line of sight 
(tt). Many groups have studied clustering of galaxies, clus- 
ters and QSOs or the cross-correlation of one with another 
using simulations and observations (e.g., Mo & White 1996, 
Mo et al. 1996, Borgani et al. 1997, Croft et al. 1999, Col- 
berg et al. 2000, Moscardini et al. , 2000, Zehavi et al. 2002, 
Hawkins et al. 2003, Croom et al. , 2005, Springel et al. , 
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2005). Springel et al. (2005) compare the galaxy 2 point cor- 
relation function from their high-resolution simulation with 
that of galaxies in the 2dF, SDSS and APM surveys and 
show that the correlation functions follow power laws for 
r <~ 20 h~ 1 Mpc and that observation and simulation agree 
very well. 

In the most recent Sloan Digital Sky Survey (SDSS) 
analysis, Zehavi et al. (2005) studied the luminosity and 
color dependence of the galaxy correlation function. It was 
found that clustering of blue galaxies increases continuously 
with luminosity, whereas bright red galaxies show strong 
clustering at large scales and faint ones mainly cluster at 
small scales. Many of the studies described fit the correlation 
function with a single power-law although this simple form 
was found be some not to be good enough to use (e.g., in 
analysis of the 2dF data by Hawkins et al. 2003, Croom 
et al. , 2005, da Angela et al. 2005). These considerations 
will affect our choice of model correlation function, described 
below. 



2.1 Theory 

Alcock & Paczyhski (1979) showed that the ratio of angu- 
lar and redshift sizes of a spherical comoving object evolve 
differently in time for different cosmologies. Detailed discus- 
sions of how the components evolve in redshift are presented 
in Ballinger et al. (1996) and Matsubara and Suto (1996). 

Da Angela et al. (2005) argued that the correlation func- 
tion (of high-z galaxies) cannot be described adequately with 
a single power law, and use a double power law when they 
fit it. We adopt this general idea when we build our model 
correlation function but instead of using two power laws, 
we use a single power-law for r ^ ro and use a correlation 
function computed using cold dark matter linear theory (the 
transfer function from Bardeen et al. 1986) for r > tq. We 
test this idea using the ACDM Hubble volume simulation 
(see §2.2 for more detailed description.) Figure Q shows the 
best-fit power law and the best-fit composite power-law and 
CDM correlation functions. It demonstrates that the latter 
describes the correlation function better, especially for large 
r where the £(r) decreases faster than a power law. 

A power law correlation function is given by: 



(1) 



r can be rewritten in terms of line of sight and transverse 
components, 



2 2.2 

r = a + 7r 

= / 2 A<9 2 + g 2 Az 2 

where / and g are given by 

f=(l + z)D A (z), 

and 

c 

9 : 



H(z) 



(2) 



(3) 



(4) 



In order to incorporate the equation of state of dark energy, 
we use (Seo & Eisenstein 2003) 
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H(z) = H \ tt m (l + z) 3 + 0.x exp 



1 + w(z) 
l + z 



dz 



(5) 
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and 



D A (z) = 



dz 



(6) 



In our study, we assume a flat universe such that fl m + 
fix = 1 and allow a redshift dependent equation of state 
with the simple form, w — u>o + W\Z (e.g., Seo & Eisenstein 
2003). We introduce a distortion parameter h (Popowski 
et al. 1998) such that 



h 



z 9 



(7) 



With respect to a particular reference cosmology, this be- 
comes 



f/g 
fo/go ' 



(8) 



To compute fo and go we use as a reference cosmology 

the same cosmology used in the Hubble volume simulation 

(fi m = 0.3, JIa = 0.7). Throughout the paper we use this 

h as the distortion parameter unless noted otherwise. Note 

that 2/£a — L j la. because A6 and Az represent the Ob- 
it ' T 9 ' 90 1 

served quantities and are invariant with cosmology. There- 
fore, h is defined in such a way that when h is greater (less) 
than 1, the correlation function is stretched in the transverse 
(line of sight) direction. 

For a given power-law correlation function used for 
r ^ ro, we adjust the amplitude of the CDM correlation 
function so that it joins the power law exactly at r = ro. An 
alternative would be to keep the amplitude of the CDM cor- 
relation function fixed but rescale r so that it joins the power 
law. This is equivalent to assuming a different value of fl/i. 
We briefly describe an analysis in §3 where this was done 
in order to test whether a better fit to the non-linear cor- 
relation function could be achieved. In our fiducial analysis, 
however we restrict ourselves to only changing the amplitude 
of the CDM £. We note that as ro determines the change in 
the amplitude, using the CDM shape does not increase the 
number of free parameters in the fit. 

In redshift space, we observe a correlation function that 
is distorted because of peculiar velocities and other noise 
sources such as redshift errors. We assume that the dis- 
tortion due to peculiar velocities can be removed using a 
measured kSZ velocity (we return to the validity of this as- 
sumption in §3.1) and we use a Gaussian function of width 
a v to model the remaining distortion. The distortions only 
occur in the line of sight direction so the effect can incor- 
porated into the correlation function using a convolution in 
■k. Using the notation of Dalton et al. (1992), the distorted 
correlation function for r ro, can be written as: 



a + 7T 



^dw (9) 



When r > ro, the CDM part of correlation function is con- 
volved in the same way: 



2tt(j v 



?CDM(r)e 2a vdw 



(10) 



where r=(a 2 + (7r-§) 2 ) 1/2 . 

Our aim is to vary h, 7, ro and o v and compare our 
model correlation function with that measured from a sim- 
ulation. This is to see if we can measure the distortion 



10° 



10' 



10' 



10" 



1Cf 



power law 
power law + cdm 
simulation 




10 



10" 



10' 

r (Mpc/h) 



Figure 1. The galaxy cluster correlation function from an octant 
of the Hubble volume simulation (points), as well as the best fit 
power law and our chosen fitting function (power law + linear 
CDM, see text.) 



parameter h. Since h is calculated with respect to the 
£l m — 0.3, SI a = 0.7 cosmology (Equation 8), we hope to 
recover h = 1, thus showing that we can recover the cosmol- 
ogy that was used to generate the ACDM Hubble volume 
simulation. 



2.2 Simulations 

We use a light-cone output of the ACDM Hubble Volume N- 
body simulation made publically available by the Virgo Con- 
sortium (Frenk et al. 2000). The output is of a model with 
flm = 0.3, S1a = 0.7, h — 0.7 and amplitude of mass fluctua- 
tions as = 0.9. The light cone data is an octant of the simu- 
lation volume in shape (and therefore covers 7r/2 steradian), 
with radius 3000 /i _1 Mpc. The redshift extends to z = 1.46. 
We use the Virgo consortium cluster catalogue generated by 
using the spherical overdensity (SO) algorithm (Lacey & 
Cole 1994) with density threshold of 200, and a minimum 
particle count per cluster of 12 (Evrard et al. 2002). The 
mass of each particle is 2.2xl0 12 h~ 1 M t r ) and there are to- 
tal of 802461 clusters with mass > 2.6 x 10 13 /i _1 M Q in the 
catalogue. We take this data set to be our fiducial "simu- 
lated survey" i.e. a survey that covers a quarter of the sky, 
down to this mass limit. 

We calculate the cluster correlation function using the 
following estimator: 



(11) 



where N cc is the number of cluster-cluster pairs and N cr is 
the number of cluster-random pairs. C is the ratio of the 
number of clusters in a random catalogue to the number of 
clusters in the cluster catalogue. Our random catalogue has 
C — 20. The random catalogues are made to have the same 
overall number density with redshift trends as the simulation 
data and the same survey boundaries. In order to do this, 
we bin the simulated clusters in redshift and use this to 
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generate a number density profile as a function of redshift. 
We then randomly sample from this number density profile 
until we reach the correct total number of random points. 

In order to measure the geometric distortion as a func- 
tion of redshift, we divide the simulated universe into 11 red- 
shift bins from z = 0.2 to z = 1.3. The space density of clus- 
ters declines at higher redshifts. For example, we have space 
densities of 3.5 x 1(T 5 , 6.0 x 10~ 6 , 3.9 x 10~ 6 clusters per co- 
moving unit volume ( /i" 3 Mpc 3 ) for z = 0.2 - 0.3, 0.7 - 0.8, 
and 1.2 - 1.3 respectively (24723, 83441 and 58767 clusters 
in the bins.) 

For each redshift bin, we compute £(<r, tv), in 2 /i _1 Mpc 
intervals. In our fitting, we keep the 494 bins in £(<r, n) which 
have r < 50 ft _1 Mpc. For \ 2 fitting, we build a covariance 
matrix using the jack-knife method with 9 subvolumes of 
equal size. We only use the diagonal elements in the covari- 
ance matrix for our analysis. This approach is intermediate 
between using Poisson errors and a full covariance matrix. 
We do not do the latter because with a small number of sub- 
samples the matrix becomes singular. Because galaxy clus- 
ters form a sparse sample, this simplification of the error 
bar calculation is reasonable (see e.g. Popowski and Wein- 
berg 1998). We check on this by computing the \ 2 P er bin 
and comparing it to unity. 
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Figure 2. The simulation correlation function §(r) for the Hub- 
ble volume octant in real and redshift space for the redshift range 
z = 0.8 — 0.9. The difference between the two in the linear regime 
agrees well with the expected Kaiser (1987) enhancement in red- 
shift space (see text.) 



2.3 Distortion in Redshift Space 

Observing clusters in redshift space means that cluster pe- 
culiar motions cause distortion in the correlation function. 
Matsubara & Suto (1996) modeled the redshift distortion of 
parallel and transverse components of the correlation func- 
tion, and Padilla & Baugh (2002) studied the redshift dis- 
tortions in the cluster correlation function using the Hubble 
volume simulation. 

Peculiar velocities reflect large-scale cluster infall into 
overdense regions and result in the correlation function be- 
ing squashed in the line of sight direction. In the linear 
regime, the flattening effect is typically described in terms 
of P = f2^ 6 /& where b is the bias parameter. The geomet- 
ric distortion and the peculiar velocity distortion both have 
roughly similar aspects, although they can be distinguished 
given large enough datasets or through their differing evo- 
lution with redshift (see e.g., Ballinger et al. 1996). 

In Figure[5]we plot the correlation function as a function 
of r in both real (£ r ) and redshift space (£ s ). We see the 
suppression of correlation function on small scales and the 
boost on linear scales as predicted by Kaiser (1987). The 
linear theory prediction is £ s /£ r = 1 + Hj^ + , where 
j3 — fi^ 6 /b, with b being the bias factor for clusters. Colberg 
et al. (2000) reported that b for these clusters is 6 = 2.25, so 
that the linear boost correction factor due to Kaiser effect 
for ACDM is 1.150. We average over the points in the linear 
regime with r > 10 /i _1 Mpc and find an average boost of 
1.146, close to the linear theory prediction. 

Although we assume that the peculiar velocites can be 
removed via kSZ surveys, we still have to deal with sys- 
tematic or measurement errors that can not be eliminated 
entirely. Unlike peculiar velocities, measurement errors will 
tend to stretch the correlation function along the line of 
sight. Padilla and Baugh (2002) showed that the cluster 
correlation function is not as subject as the galaxy corre- 



lation function to small scale virialized motions ( "the finger 
of God" effect), due to the fact that superclusters are not 
virialized systems. 

In Figure|3]we examine £((T,7r) measured from the sim- 
ulation. We show how the correlation function is distorted in 
redshift space when there are random errors (two right pan- 
els) and when there are peculiar velocities (two bottom pan- 
els). For the random errors, we add Gaussian redshift errors 
with a — 400km/s randomly to each cluster. The top left 
panel shows the case when there is no redshift distortion at 
all. The correlation function is isotropic, as expected. With 
random errors only, it is stretched along the line of sight 
(top right), and undergoes further distortion when peculiar 
velocities are added (bottom right). The overall squashing 
along the line of sight is better seen in the bottom left panel 
where peculiar velocities are only source of distortion (no 
random errors added). This agrees well with the results of 
Padilla & Baugh (2002) who find that there is no stretching 
of the correlation function due to random peculiar veloci- 
ties in redshift space, but that coherent motions of clusters 
flatten the correlation function along the line of sight. 

Nagai et al. (2003) studied the real and estimated pecu- 
liar velocities and their discrepancies with a hydrodynami- 
cal simulation. We describe their work in more detail in §2.5 
below. We take their results to be representative of the min- 
imum noise level which can be achieved in measurements of 
the cluster velocity. 

2.4 SZ effect 

The SZ effect comprises the thermal and kinetic effects. The 
thermal SZ effect is a distortion of the CMB photon energy 
spectrum due to inverse Compton scattering in the cluster. 
When the CMB photons are scattered by free electrons in 
the hot intracluster gas, they exchange energy and the tem- 
pature decreases (increases) at low (high) frequencies, result- 
ing in a spectral distortion skewed toward higher frequencies. 
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Figure 3. Measured correlation functions from the Hubble Vol- 
ume simulation for the redshift range z = 0.8 — 0.9. Left (right) 
panels: without (with) Gaussian measurement errors of 400 km/s. 
Top (bottom) panels: without (with) peculiar velocities. The cor- 
relation function is almost isotropic when there is no redshift 
distortion (top left). 



The distortion due to the thermal effect (non-relativistic) is 
described by: 

AT tsz ( e* + l 



TcMB 

where x 
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In addition, when the cluster is moving with respect 
to the CMB rest frame, the signal will be Doppler shifted, 
with an amplitude proportional to the peculiar velocity of 
the cluster. This is the kinetic SZ effect. The magnitude of 
the effect is given by 

AT kS z 



TcMB 



I Vpec \ 
\ C J 



(14) 



where r e is the optical depth and v pec is the peculiar veloc- 
ity component along the line of sight. The clusters moving 
toward (away from) the observer have a negative (positive) 
sign for the temperature decrement. For detailed reviews 
of the SZ effects, see e.g., Sunyaev & Zel'dovich (1980), 
Rephaeli (1995), Birkinshaw (1999), Carlstrom et al. (2002). 

As Holder (2004) emphasized, only 1 % of CMB photons 
interact with free electrons in the hot cluster gas and the en- 
ergy exchange per scattering is another 1 %, which makes the 
total thermal effect on the order of 10 -4 . The kinetic effect is 
one magnitude smaller than this (10 ). Although both ef- 
fects are very small, they have different spectral signatures. 
For the thermal SZ, the intensity decreases (increases) at 
low (high) frequencies, vanishing at ~ 220 GHz, whereas at 



the same frequency the kinetic effect is at its peak and it has 
constant sign (positive/negative) throughout all frequencies. 
With a multiwavelength SZ survey, we expect to be able to 
measure both the thermal and kinetic effects, from which we 
can determine three physical properties of each cluster (t, 
v pec and T e ). 

The cluster temperature can be obtained from X-ray 
observations or from thermal SZ measurements when the rel- 
ativistic corrections are included (Pointecouteau et al. 1998, 
Hansen et al. 2002). Sehgal et al. (2005) pointed out how- 
ever that the SZ measurements are not sufficient to break 
the degeneracy between the 3 cluster parameters (r, v pec 
and T e ) and showed that an independent measurement of 
cluster temperature would greatly help. However, X-ray ob- 
servations may not guarantee accurate information on the 
cluster temperature either. One difficulty stems from the 
fact the X-ray emission signal decreases with distance due 
to the inverse square law, making measurements at higher z 
much harder (unlike the SZ effects). An additional problem 
arises because in inhomogeous cluster gas, the temperature 
inferred from X-ray emission is affected by the dumpiness 
of the intra cluster gas while it is irrelevant in thermal SZ 
measurements (e.g., Hansen 2004). 

A recent study of systematic effects was carried out by 
by Diaferio etal (2005), who used hydrodynamic cosmologi- 
cal simulations to create large samples of simulated clusters. 
It was found that it is crucial to use the electron weighted 
cluster temperature, T e to recover the peculiar velocity from 
the kSZ effect. Using the X-ray emission-weighted tempera- 
ture, Tx can overestimate the peculiar velocity by 20 — 50%. 
Spatially resolved nearby clusters can be used to measure T e 
in the center where it is comparable to Tx- However, spatial 
modelling of the X-ray emission is still needed to seperate 
the SZ effects, resulting in potential overestimate of the pe- 
culiar velocity by 10 — 20%. 

The accuracy of v P e C measurements depends on how 
well one can reduce the systematic errors. Many studies 
(e.g. Knox et al. 2004) have been carried out testing how 
well one can hope to extract v pe c via the kSZ effect, in- 
cluding the effects of known systematic errors and noise 
sources, e.g., interstellar dust emission, infra-red (IR) galax- 
ies, radio sources, imperfect bulk velocity, etc. The primary 
CMB fluctuations which also impede accurate measurement 
of the kSZ because they have the same spectral behavior 
(see Haehnelt & Tegmark 1996, Aghanim et al. 2001, Forgi 
& Aghanim 2004, Aghanim et al. 2005). Some studies have 
estimated the potential measurement error on v pe c to be as 
low as ~ lOOfcms (Nagai et al. 2004, Holder 2004, Sehgal 
et al. 2005), while Benson et al. (2003) have set an obser- 
vational upper limit on cluster velocities with errors of the 
order of ~ 1500fcms. Our study is based on the former esti- 
mates, however we have also tried ~ 2000fcms as a measure- 
ment error. In this case we found that the error bars were 
too large to be useful for constraining cosmological param- 
eters and we could not recover the distortion parameter (h) 
sensible enough to continue for further study. If the mea- 
surement errors are too big, we can not tell wheather the 
distortion is due to A or v pec - This may be the most com- 
mon but challengin problem in AP test along with the fact 
that we need so many accurate measurements of the clusters 
velocities. 
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2.5 kSZ Errors from Hydrodynamic Simulations 

The measurement of peculiar velocity from the kSZ effect 
in observational data will be affected by many sources of 
noise. In our model fitting, we use a Gaussian function to 
parametrize the statistical scatter on measurement results. 
In creating the simulated surveys from the Hubble volume 
simulation we will include the option of using a more realistic 
noise distribution. In this section, we use the error distribu- 
tions from the hydrodynamic simulation work of Nagai et al. 
(2003). 

These authors use a high-resolution simulation of a 
galaxy cluster by Kravtsov (2002) to generate detailed kSZ 
maps in a ACDM cosmology with Q m = 0.3, SIa = 0.7, h = 
0.7. This is an iV-body+gasdynamics simulations, but does 
not include the effects of magnetic fields, gas cooling, stellar 
feedback or thermal conduction. Nagai et al. use their sim- 
ulation to create simulated kSZ maps of the cluster. They 
then analyse these simulated maps in the same way one 
would with real observational data, and use this to com- 
pute an kSZ-inferred cluster peculiar velocity. This value is 
then compared to the known peculiar velocity of the clus- 
ter in the simulation taken from the particle data. The idea 
behind the paper is to see how well the two values match 
up, what the scatter and systematic differences between the 
two are, and how the measurement can be optimized, For 
example, for the simulated kSZ maps, Nagai et al. try cal- 
culating the density-weighted velocities within certain radii. 
They find that different aperture sizes do lead to different 
results for the scatter between the true and inferred veloc- 
ity. This is because clusters are not perfectly spherical and 
have internal structures and motions. If the aperture size is 
too small it may not reflect the bulk motion of entire cluster 
properly. We use the data calculated using the virial radius, 
which they find is the most optimal case. 

A set of estimated and observed values of kSZ velocities 
was kindly provided by Daisuke Nagai. The simulated obser- 
vations were generated for three orthogonal projection an- 
gles and at nine different epochs (a = 0.60, 0.65, 0.95, 1.0 
where a = 1/(1 + z).) We will assume that these orthogonal 
projections represent independent measurements. The pecu- 
liar velocity measurement errors are the difference between 
the estimated and observed radial velocities. The distribu- 
tion of Vr Bt — v° 3 that we use is shown in Figure 6 of Nagai 
et al. . The standard deviation of these measurement errors 
is a — 50 km s . 

For each of the clusters in our Hubble volume simula- 
tion catalogue, we randomly pick a value from this set of 
adopted measurement errors, at the closest appropriate red- 
shift. We then add it to the cluster's line of sight Hubble 
velocity and compute the new redshift. We next convert the 
cluster coordinates into comoving coordinates assuming an 
LCDM cosmology and compute the correlation function as 
normal. As a result, we have a cluster correlation function 
distorted by measurement errors, one that would be seen 
from a survey if the only source of error was due to dif- 
ferences between the projected and actual cluster velocity. 
This represents the best case that could be obtained from a 
survey. In §4, however, we will use larger and more realistic 
Gaussian errors when we simulate ACT or SPT-like surveys. 

As mentioned above in reference to the work of Diaferio 
et al. 2005 and others, systematic biases in measured veloc- 
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Figure 4. Contours of Ax 2 for fits of the parameters h,ro,7 and 
<j v for z=0.8-0.9. We fit the 2 point correlation function from the 
Hubble volume simulation to the power law + CDM correlation 
function I - Eanl9b 1101 . The measurement errors from the hydro- 
dynamic simulation of Nagai et al. (2003) have been added. A\ 2 
in 4 parameter space has been marginalized for the two parame- 
ters of interest in each panel. Contour lines represent 1, 2 and 3 
a confidence levels on each of 2 parameters. 



ities, such as over or underestimates are also expected to 
occur in real measurements. Although in our fiducial mod- 
elling, we use the errors from bulk flow modelling, we will 
also estimate the effects of systematic biases on cosmic ge- 
ometry. 



3 METHOD 

In order to constrain cosmic geometry from the simulated 
cluster catalogue, we fit our model correlation function 
(§2.1) to it, varying four free parameters, h, ro, 7, and a v . 
The parameter ro determines the overall amplitude of the 
correlation function, and 7 its slope. The line of sight veloc- 
ity distribution added to model measurement errors is given 
by a v . The geometric distortion, stretch in the transverse 
direction is quantified by h. Our purpose is to measure the 
distortion parameter h and its uncertainty as a function of 
redshift by marginalizing over the other three parameters. 
We then use h(z) to find the corresponding cosmology. The 
detailed procedure is explained in the next section. 

3.1 Model fitting 

The model correlation function will be isotropic when there 
is no redshift distortion due to peculiar velocities or mea- 
surement errors. However, the measured correlation func- 
tion will not be. The two more parameters in addition to 
ro, 7 which we have introduced will parametrize the isotropy 
of £(cr, 7r). This is done by fitting the measured correlation 
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Figure 5. A\ 2 for fits of the geometric distortion parameter 
h, and the two correlation function parameters r and 7 after 
marginalization over the other three parameters in each case (see 
§3.1). 



function with the distortion parameter h. In principle, this 
should yield h = 1. However, sampling variations will yield 
a measured correlation that is not perfectly isotropic. The 
statistical uncertainty on h measured from a particular cat- 
alogue will depend largely on these sampling variations. 

The fitting is done for z = 0.2 ~ 1.4 and r = 1 ~ 
50 /i -1 Mpc and the procedure is as follows: 

(i) Add measurement errors to the cluster redshifts in the 
simulated catalog. 

(ii) Convert angular positions and redshifts of clusters to 
comoving coordinates assuming an LCDM cosmology. 

(iii) Calculate the correlation function. 

(iv) Compute a model correlation function parametrised 
by values of h, ro, 7, a v 

(v) Calculate y 2 for the fit using the covariance matrix. 

(vi) Marginalize the error in the parameter of interest. 

(vii) Find the best fit parameters that yield the minimum 

x 2 . 

Since we calculate \ 2 with 4 parameters, we need to 
marginalize it for each parameter. For example, to marginal- 
ize % 2 m h we rewrite x 2 m terms of the likelihood, 

-21n£ 



X 
then 

C — exp" 



(15) 



(16) 



with this, calculate the average likelihood over ro, 7 and o v 

. , x 2 

£ = \ exp t Ar A7Ao\, (17) 

ro,7>o"„ 

where Ar=0.02 /i _1 Mpc, A7=0.005 and Acr„=2.5 km s" 1 . 
Now we calculate x 2 f° r h from jC 



2 

Xh 



-21n£. 



(18) 



In order to demonstrate our fitting procedure, we have 
applied it to one redshift bin in a simulated survey. This 
survey is the fiducial one, including all the clusters in the 
Hubble volume octant. The particular redshift interval we 
choose is z = 0.8 — 0.9 and includes ~ 80000 clusters. In 
this example, we have applied the hydrodynamic simulation 
velocity errors only to the cluster redshifts, so that we expect 
to recover a low value of the velocity dispersion parameter 
a. 
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Figure 6. Measured (left) and best fit (right) correlation func- 
tions for the Hubble volume octant between redshifts 2 = 0.8 and 
2 = 0.9 (this is ~ 80000 clusters). Measurement errors from the 
hydrodynamic simulation have been added. The effect of these is 
not as visible as in Figure U] where much larger errors are added 
(400 km/s). For the bottom panels, peculiar velocities are also 
included in addition to the measurement errors. 



in two-parameter space, with separate plots for each pair of 
parameters. In each of the panels, the curves represent the 
Ax 2 values after marginalization over the other two param- 
eters. The one dimensional distributions, from A^ 2 values 
marginalized over the 3 other parameters (equation 1171 arc 
shown in Fig. We find la statistical errors on h, ro and 
7 of 0.5%, 0.4%, and 0.3% respectively. The important pa- 
rameter for our purposes is h, for which we find a best fit 
value of 0.987, which is 1.3% (~ 2.5cr) from the expected 
h — 1. We attribute this percent level bias to the fact that 
the model correlation function we fit with is not a perfect 
match to the simulation correlation function. This type of 
systematic bias is likely to be difficult to circumvent. We re- 
turn to this in the discussion (§5.2). The chi squared value 
for the best fit in this case is 578, for 494 bins, so while the 
fit is good, there is room for slight improvement. 

We examine the £(<r, n) plots from the simulation in 
Fig.|S| alongside the best fitting model £. In the top panel of 
Fig. |B]we have included only the measurement errors from 
hydrodynamic simulations, but assumed that the peculiar 
velocities have been subtracted. Although the velocity error 
distribtion is not Gaussian, we have fitted it in our convo- 
lution using a Gaussian function (equations l9l fc [TUt . which 
has a best fitting a ~ 50km/ s (see figBJ. By eye, it is diffi- 
cult to see effect of these velocity errors making the £(<r, n) 
anisotropic, although this effect sucessfully recovered in the 
fitting. 

In order to judge the effect of peculiar velocities on a 
straight fit to the cosmic geometry, we have also tried fitting 
our 4 parameters without subtracting the peculiar velocities. 
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By doing this, we will get an idea of how the peculiar veloc- 
ities can mimic the effect of geometric distortions. This was 
addressed in the context of linear theory by Ballinger et al. 
(1996), who found the relationship governing the effective 
value of h is h ~ 1 + 2/3. Here f3 — Q.®f /b, with b being 
the bias parameter. In the present case, b (the ratio of the 
cluster £ to the matter £) is ~ 2.5, so that /3 ~ 0.2 — 0.4 over 
redshifts from z = — 1. Based on this, we expect h to be 
~ 1.4 — 1.8 when peculiar velocities are not removed. 

This simulation test was carried out, and £(<r, n) is 
shown in the bottom panels of Fig. [5] alongside the best 
fitting model which includes geometric distortions only. We 
can see that the peculiar velocities cause a very strong dis- 
tortion, and that the form of the distortion is not obviously 
different to the eye when compared to the geometric one, 
a point which has been made by many authors including 
Ballinger et al. (1996) and Matsubara and Suto (1996). 
When we fit the h(z) that results (Fig. |7J we find values 
~ 1.6-1.8, as expected. The x 2 per bin is 3.2 (1616 for 494 
bins), which is substantially worse than the fit to true geo- 
metric distortions. This difference in the fit could be useful if 
accurate estimates of the peculiar velocities are not removed. 

For example, a systematic overestimate of the velocity 
will cause the distortion h(z) to be underestimated. In the 
worst case scenario detailed by Diaferio et al. 2005, using a 
naive measurement based on clusters unresolved in X-rays 
would result in an overestimate of the peculiar velocity of 
20 — 50%. Subtracting these cluster velocities to recover an 
estimate of the real space correlation function would yield 
h(z) values ~ 0.1 — 0.35 too low. This is at the level which 
would make an Einstein De Sitter model look like a concor- 
dance A model, and so this would not result in an acceptable 
test of cosmic geometry. Even with a spatially resolved mea- 
surement of the temperature, from Diaferio et al. 2005 we 
would expect a bias in h(z) of ~ 5 — 10%, which would 
affect a measurement of Qa at the ~ 0.1 — 0.2 level. We 
note that Seghal et al. (2005) find somewhat more opti- 
mistic conclusions from a simulation study of constraining 
cluster parameters from SZ observations, finding that ve- 
locity errors could be biased by as low as 15 — 40 km s^ 1 
if X-ray weighted temperature measurements are used. This 
said, there are also some improvements which could be made 
to reduce the bias, such as using the integrated SZ flux to- 
gether with scaling relations (Benson et al. 2004), and/or 
making use of the z dependence of the distortion more effi- 
ciently. As we have seen above, it is also possible that the 
relatively poor goodness of fit when peculiar velocity dis- 
tortions remain will help diagnose this without a theoretical 
model for the peculiar velocities being necessary. 

With presently planned surveys such as ACT and SPT, 
these biases will be comparable to or smaller than the statis- 
tical errors, however they will represent an severe obstacle 
to any attempt to use the cluster real space AP test as a 
precise probe of cosmology. 

In Fig. |7| we also show h vs. z from our fitting to the 
simulation after subtracting peculiar velocities, but with the 
errors from hydrodynamic simulations. We can see that the 
expected result for LCDM, ft = 1 is recovered well, with 
error bars of ~ 1%. The trend of h(z) expected in some other 
example cosmologies is also shown as lines. The differences 
in h(z) between an Q. m = 1 model and LCDM are of the 
order of 10s of percent. Models with slightly different dark 
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Figure 7. Distortion parameter h vs. z for different cosmologies 
with repect to an Q m = 0.3, Qa = 0.7 cosmology. We show the- 
oretical results for a few example cosmologies as lines (Q m = 0.3 
unless noted otherwise). For the simulations, we show results as 
circles (with 1 a error bars) for the case with peculiar veloci- 
ties corrected when random errors from hydrodynamic simulation 
are added. We also plot as triangles results for the simulations 
where much larger velocity errors were added (a Gaussian a of 
400 km s" 1 ) Finally, we also show results for the case where the 
peculiar velocities have been left in, as crosses, which yields an 
extremely distorted correlation function (h very different from 1.) 



energy equations of state from LCDM fall even closer. We 
explore these below. 

We have also carried out simulation tests using much 
larger peculiar velocity errors. For example, we have as- 
signed random Gaussian distributed errors of 400 km s' 1 to 
the clusters before carrying out our fitting procedure (Fig- 
ure 3). In this case, we find values of h(z) which are approxi- 
mately 3% too high, when averaged over the 10 redshift bins. 
This corresponds to a quite large bias, and means that it 
seems to be difficult to recover the correct cosmic geometry 
with such large velocity errors. With errors ~ 200 km s _1 , 
we find that the fitting technique works much better. Er- 
rors of this magnitude are to be expected given good ob- 
servational data (see e.g., Diaferio et al. 2005, Seghal et al. 
2005). We have not so far mentioned the errors on the clus- 
ter redshifts, but of course an additional complication is that 
they must also be measured to a comparably high level of 
precision. 



3.2 Varying the dark energy equation of state 

We adopt a simple equation of state, w = wo + Wiz, where 
the pressure, p = wp. In order to find the best fit values 
for wo and wi, we vary wo, ttfi and Qm on a grid, and 
compute the expected h(z) for each set of parameters. We 
then calculate Ax 2 by fitting to the simulation results in 
each redshift bin. \ 2 ls given by: 



* 2 = E 



(hi - Kf 



(19) 
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Figure 8. Fitting cosmological equation of state parameters (see 
§3.2) to the fiducial Hubble Volume simulation catalogue. Best 
fit at wo = —1.08,toi = —0.15. We show 1,2 and 3cr contours. 
The top panel is without a prior on Q m and has a minimum 
at O m = 0.32. The bottom is with a prior of Q m = 0.3, with 
uncertainty on Q m of a = 0.05. For visual purposes, a Gaussian 
filter with a = 0.8 grid cells was used to smooth before plotting. 
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Figure 9. Distortion parameter h(z) measured from the fiducial 
Hubble Volume simulations (points) together with the theoretical 
curve for the best fit wo,wi and H m (which gives the minimum 
X 2 in Figure 8.) The line for the correct model ACDM is h(z) = 1. 
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where hi represents the best-fit distortion parameter h in 
redshift bin i, h'f is the distortion parameter with given wo, 
Wi and CIm, and the sum is over n bins in redshift. In Eg n il 91 
the cr's have been obtained from the ID marginalized results 
for h(z) from the simulations (e.g. Figure 5). In the few cases 
where in the higher dimensional space before marginaliza- 
tion, the la contour for h is out of the range, we are conser- 
vative and set the \ 2 to a very large number so that during 
the marginalization the contribution to the likelihood from 
that point is negligible. 

Once we have a 3d grid of A\ 2 values for Wo,Wi and Qm 
we marginalize over Q m . We do this either with no external 
prior, or else after imposing a Gaussian prior on Q m (C = 

(fi,„-0.3) 2 

C exp c os ). We expect to recover w$ = — 1, w\ = 0, as 
this is the cosmology used in the Hubble volume simulation. 

From our fitting procedure we find wo = — l-OS^ol! 
and wi = — 0.15g|i7 with no prior imposed on Q m (la error 
bars). Adding the prior makes little difference, with the best 
fit being w = -1.09^26 and wi = O.lg;^. In Figure 8 we 



Figure 10. Constraints on wq and Q m for the fiducial Hubble 
Volume simulation catalogue case when we set w\ = 0. We show 
1, 2 and 3a contours. The cases with and without a prior on fim 
give almost the same result. Best fit at wo = —1.09 and fi ra = 
0.31. 



show the Ax 2 contours for these two cases. Although we 
have seen from Figurc|H]that the distortion parameter h can 
be recovered at the percent level, this still translates into 
relatively loose constraints on W\ , regardless of the prior on 

In Figure [5] we show the shape of the theoretical h(z) 
curve for this best fitting cosmology, compared to the mea- 
sured values from the simulation. We can see that the best fit 
curve follows the dip in h(z ) values at high z, which account 
for the difference from the true cosmology. 

In Figure Hill we set w\ = 0, in order to assess the type 
of constraints possible in this case. We find a minimum x 2 at 
Q m — 0.3l[5;oo2, and wo = — 1.09o:ol?- Imposing the previous 
prior on Q m makes no difference to this result. 
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Figure 11. Theoretical distortion parameter h(z) for tCDM and 
ACDM (lines). Data points with circles (crosses) represent fitted 
distortion parameters from the simulation with 1 <r error bars for 
rCDM (ACDM). 



3.3 Test on the tCDM simulation 

We now test our method on a different simulated universe. 
In the previous case, have taken data sets from ACDM sim- 
ulations, and then assumed the ACDM geometry when an- 
alyzing them. Because of this, we expected to recover a dis- 
tortion parameter h(z)=l. Now instead, we use a different 
simulation, but still assume the ACDM geometry to analyze 
it. In this case, we expect to find h(z) 7^ 1. 

For the simulated universe, we use the Hubble volume 
simulation of the rCDM model (see Frenk et al 2000). This 
is a cosmology with Q m = 1,Qa = 0. The linear matter 
power spectrum has a similar shape to that of LCDM. We 
use all the clusters in an octant as before, except this time 
the maximum cluster redshift is z max = 1.25. 

We first compute the angular positions and redshifts of 
the clusters, as they would be seen by an observer. We then 
use the ACDM relations (Equations 5 and 6) to convert 
these into comoving coordinates, assuming an observer at 
the origin. We compute the correlation function from the 
simulation as before, and fit the distortion parameter h(z). 

This test is is analogous to the real situation, where we 
would have observational data from an unknown cosmology 
and use an assumed cosmology during the fitting procedure. 
The unknown cosmology here is an Einstein-de Sitter uni- 
verse with Q, n — 1, and we find the distortion parameter 
with respect to a cosmology with Q. m = 0.3, Qa = 0.7. 

Using Equation 8 we compute the distortion parame- 
ter which we expect to recover in this case. This is shown 
as a dotted line in Figure [TT1 When we carry out the test, 
we find results which are intermediate between this line and 
the h(z) = 1 LCDM line at all redshifts. This is an indica- 
tion that in a realistic situation such as our test, the distor- 
tion parameter cannot be inferred simply when the assumed 
model (here LCDM) is very different from the actual Uni- 
verse (in this case rCDM) . Instead, if the distortion param- 
eter is found to be very different from h — 1, an iterative ap- 



2 

Figure 12. Distortion parameter h(z) with 1 cr error bars for the 
simulated SPT cluster survey. 3 example theory curves are shown 
as lines. 



proach should be used. In this case, the assumed cosmology 
on the second try would be one intermediate between LCDM 
and tCDM. By using this method, we would approach the 
true cosmology, albeit more slowly. In Figure HTI we show the 
recovery of the Einstein-de Sitter geometry as points with 
error bars, using the EDS model as an assumed cosmology. 
This is representative of what should occur with the full it- 
erative procedure, although this has not been fully tested in 
this case. We have rescaled the result for h(z) by the ratio of 
the analytical predictions for ACDM and rCDM so that we 
can show the results on the same plot as the ACDM results. 
Although the fit is not quite as good as with the ACDM 
test, the results show that by moving towards an estimate 
of an assumed cosmology we can expect good results. 



4 RESULTS FROM SIMULATED 
OBSERVATIONAL CATALOGUES 

4.1 The South Pole Telescope 

The SPT survey will cover 4000 square degrees in 5 fre- 
quency bands and is expected to observe ~20,000 clusters 
with masses greater than 2xlO 14 M0. For our SPT (and 
ACT, see below) simulations, we use a Gaussian measure- 
ment error of 100 km/ 's, based on published estimates of the 
best likely error on the kSZ velocity measurements (Nagai 
et al. 2004, Holder 2004, Sehgal et al. 2005). The standard 
deviation of errors from the hydrodynamic simulation which 
we used before in the Hubble volume case was ~ 50km/ s, so 
the present value gives slightly bigger velocity distortions. 

The simulation output we use to make our mock cata- 
logue is an octant shape and covers almost the same area (we 
again use the ACDM simulation), so to simulate the SPT 
survey we set a mass threshold so that the number of clusters 
that we observe is around 20,000. This means that the total 
number of clusters in the mock survey is only 1/40 of that in 
the underlying simulation. Since we have fewer clusters, we 
increase the size of each redshift bin and include clusters up 
to z = 1 only. Our covariance matrix is constructed in the 
same way as in §3.1. When building the covariance matrix, 
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Figure 13. Distortion parameter h(z) with 1 cr error bars for the 
simulated ACT cluster survey. 3 example theory curves are shown 
as lines. 
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Figure 14. Constraints on Qa for mock ACT and mock SPT 
surveys, assuming that w\ = and Q m + Qa = 1. 




we make sure that a, tv bins are large enough that there are 
cluster pairs in each bin, to avoid singular error bars. 

The results are shown in Figure 1121 As expected, we 
have large error bars because of the of the smaller number 
of clusters than in the fully sampled simulation. The points 
also lie somewhat above the h(z) = 1 curve, but the results 
still give a reasonable estimate of cosmic geometry without 
a large bias. Due to the large error bars, it is not possible to 
constrain wq-wi within reasonable error bounds. However, 
by assuming that the dark energy density does not change 
with redshift, we can still find useful constraints on Q. a, and 
wo, for example. In Figure tTll we show the constraints on S7a 
which result when a flat cosmology is assumed and wi = 0. 
We find 1 sigma errors on Q,a of lg oli ano - a central value 
within 1 cr of the correct one. The constraints on wo (Figure 
1151 are weaker, with wq < —0.76 at 1 a. 



4.2 The Atacama Cosmology Telescope 

The ACT will cover only 200 square degree of the sky, but 
at higher sensitivity than the SPT. It is expected to observe 
~ 1000 clusters with masses greater than 2 x 1O 14 M0. With 
this mass threshold we find fewer cluster in the Hubble Vol- 
ume simulation in a comparable area, so we use 1O 14 M as 
our threshold. The x 2 analysis has been done in the same 
way as before with a Gaussian measurement error of 100 
ferns, except that the covariance matrix has been built us- 
ing 6 separate volumes of 200 square degrees taken from the 
Hubble Volume octant. We then take one subvolume as our 
simulated survey and carry out our % 2 analysis using that. 

Since we now have even fewer clusters than in the mock 
SPT survey, the error bars are substantially larger (Figure 
1131) . It is not possible to put any useful constraints in the 
wq-wi plane without other priors. Again we we try to see if 
it is feasible to constrain Qa setting Wi = (Figure 1141 and 
wo (Figure 15). We find 1 a errors on Qa ~ 2 times as large 
as for the mock SPT survey and a lcr limit on Wo < —0.45. 

From the analysis of ACT mock surveys we can see that 
the relatively small cluster sample and the need to split it 
up into redshift bins will make it difficult to measure the 




Figure 15. Constraints on wo for mock ACT and mock SPT 
surveys, assuming that w\ = and Q m + Qx — !• 



correlation function accurately enough to make this type of 
measurement. Nevertheless, some constraints will be possi- 
ble, and the error bars of e.g. ~ 0.2 — 0.3 on Qa are still 
useful. 



5 SUMMARY AND DISCUSSION 
5.1 Summary 

Using a combination of the Hubble Volume simulation 
(Frenk et al. 1999) and the hydrodynamic simulation re- 
sults of Nagai et al. (2003), we have tested using the pecu- 
liar velocity corrected galaxy cluster correlation function to 
measure cosmic geometry. Our findings can be summarised 
as follows: 

(1) The galaxy cluster correlation function, corrected using 
kSZ velocities can act as an intrinsically isotropic object for 
the AP test. 

(2) The geometric distortion of the correlation function 
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needs to be measured to % level accuracy in several red- 
shift bins to constrain the equation of state parameter wo 
to better than 10%. 

(3) Our simulation tests assuming errors from hydrodynamic 
modelling show that this can be done with a survey covering 
l/8th of the sky, with clusters out to z ~ 1. 

(4) The distortion from peculiar velocities will be very diffi- 
cult to correct. Systematic over or underestimates of the pe- 
culiar velocities from kSZ measurements will be more prob- 
lematic than statistical errors. 

(5) Of surveys currently underway or planned, the SPT 
seems to be the most promising, and could yield constraints 
on Qa accurate to 10 — 20%. 

5.2 Discussion 

Of the many recent papers which have been written on us- 
ing the AP test (e.g. Mc Donald et al. 2003), all take the 
approach of modelling the peculiar velocity component of 
the correlation function distortion. The fact that for galaxy 
clusters the peculiar velocity can be instead measured and 
subtracted offers the potential of a useful complementary 
method. For example, in order to model peculiar velocities 
well, the relationship between the clustering of the objects 
and the underlying mass field needs to be known to high pre- 
cision. This is not necessary in the case where the velocities 
are subtracted. The latter method does however bring with 
it several new difficulties, including the necessarily sparse 
nature of galaxy cluster samples and the extremely chal- 
lenging problem of measuring the peculiar velocities of thou- 
sands of clusters to high accuracy, when to date this has not 
yet been done for a single cluster. 

One of the ingredients in our measurement of the distor- 
tion is our model for the correlation function as a function 
of separation r. We have found that using a fitting function 
which does not model the correlation function well can lead 
to biases in the measured distortion. For example, before 
adopting our fiducial fitting function, a power law combined 
with a CDM linear theory £, we tried a using single power 
law as the fitting function. This resulted in a systematic 10% 
bias in the value of h(z). We also tried using a single power 
law to fit £ to regions r — 1 ~ 40 /i _1 Mpc (following da 
Angela et al. 2005). In this case the fitted h in each redshift 
bin changes slightly but the overall fitting of cosmological 
parameters does not improve very much. Using a more so- 
phisticated fitting function than a power law therefore seems 
to be crucial. 

The spectral distortion in the SZ effects is determined 
by r, v pec and T e . These three parameters could be disentan- 
gled using multiwavelength observations. Sehgal et al. (2005) 
have shown how this will be possible in light of current sur- 
vey sensitivities and suggest an independent temperature 
measurement from X-ray observations in order to break de- 
generacies. However, as mentioned in §2.5, the temperature 
obtained by X-ray spectroscopy for unresolved clusters is 
not the electron- weighted temperature which is required to 
separate three parameters (r, v pec , T e ). Mathiesen & Evrard 
(2001) estimated that the discrepancy can be as much as 
IkeV. However, Sehgal et al. found that this temperature 
discrepancy will lead to a velocity bias of 15-40 km/s. They 
also claim that for an ACT-like survey, a 2 keV temperature 
error can still only lead to peculiar velocities with errors 



less than 100 km/s. We have found that the statistical mea- 
surement errors of this size can be modelled well, but that 
any systematic bias in the velocities will give biased results 
for cosmic geometry. This issue will need to be explored in 
more detail once large samples of clusters with kSZ mea- 
surements are available. Indeed, if the bias is too large, then 
it is possible that by assuming a model for cosmic geometry 
the distortion can be used instead to calibrate the peculiar 
velocities. 

There many other planned surveys than ACT and SPT 
that are designed to be sensitive to SZ signals from clusters. 
For example, AMiBA 4 , AMI 5 ( ~ 100 clusters), APEX 6 
(~ 1000 clusters) and Planck 7 (~ 10000 clusters). The 
Planck mission promises the largest survey area (the whole 
sky), compared to SPT for example which will observe 
20,000 clusters over 4000 square degrees (about 1/8 of the 
sky). Planck will however have much larger CMB tempera- 
ture measurement uncertainties, and so will only observe a 
similar number of clusters. They will therefore be distributed 
even more scarcely, making correlation function measure- 
ments difficult. We have tried sampling our simulation down 
to the space density of Planck clusters, but find that the 
Poisson errors on £ make it unfeasible to extract useful con- 
straints on the geometric distortion. 

In conclusion, we have shown that by subtracting kSZ 
measured velocities, it is possible to carry out a test of cos- 
mic geometry with an intrinsically isotropic object, as envis- 
aged in the original paper by Alock and Paczynski (1979). 
Using simulations we have seen that given small enough sys- 
tematic and statistical errors, the correct geometry can be 
recovered. However, huge samples of galaxy clusters are re- 
quired with precisely measured velocities. For example, to 
measure wq with 3% statistical errors we have seen that 
800,000 clusters are necessary. As the kSZ signals are so 
small, observing the kSZ effects from this many clusters is 
not possible with currently planned surveys. However more 
advanced surveys are planned to detect more clusters in the 
future and the test we have proposed here will offer an ad- 
ditional way to use such data to constrain dark energy. 
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